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I. INTRODUCTION 

Quantum walks in several topologies |l| are being studied as potential sources for new 
quantum algorithms. Recently, quantum search algorithms based on different versions of 
the quantum walk, have been proposed 

as 

. These algorithms take advantage of quantum 



in the context of quantum walks. The first studies jj, |5( where numerical and considered 
walkers driven by two coins which where maximally entangled by their initial condition. 
More recently, the coin-position entanglement induced by the the evolution operator of a 
quantum walk on the line was investigated numerically. For pure states, entanglement can 
be quatified using the von Neumann entropy of the reduced density operator (entropy of 
entanglement). Using this measure, it was suggested [6] that for all coin initial states of 
a Hadamard walk the entanglement would have a long-time value close to 0.872. In this 
work, we use the Fourier representation of the Hadamard walk on a line to investigate the 
dependence of asymptotic entanglement with the initial conditions. This powerful approach 
enables us to derive the condition that must be satisfied by the initial state in order to 
obtain a limiting value of 0.872. However, for localized initial states with generic coins the 
asymptotic entanglement changes smoothly between ~ 0.736 and 1 (maximum entangle- 
ment). Furthermore, if non-local initial conditions are considered, all entanglement levels 
are accessible in the long-time limit. We consider in detail the case of non-local initial 
conditions restricted to a certain position subspace. 

This work is organized as follows. In Section [TTJ we briefly review the discrete-time 
quantum walk on the line and define the entropy of entanglement. To illustrate the general 
method, the asymptotic value for the entanglement of a particular localized initial condition 
is obtained analytically. An expression for the dependence of the asymtptotic entanglement 
for generic localized initial conditions is derived in the Appendix. Section II I II is devoted 
to the calculation of the asymptotic entanglement induced by the evolution operator of a 
Hadamard walk for a particular class of non-local initial conditions. Finally, in Section IIVI 
we summarize our conclusions and discuss future developments. 
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II. DISCRETE-TIME QUANTUM WALK ON THE LINE 



The discrete-time quantum walk can be thought as a quantum analog of the classical 
random walk where the classical coin flipping is replaced by a Hadamard operation in an 
abstract two-state quantum space (the coin space). A step of the quantum walker consists 
of a conditional traslation on the line. Of course, if the quantum coin is measured before 
ta^g a step, the ciaasrca! w* is recovered fl ^ Tire Hilbert space H = H P * H c , is 
composed of two parts: a spatial subspace, TCp, spanned by the orthonormal set {\x)} where 
the integers x = 0, ±1, ±2 . . . are associated to discrete positions on the line and a single- 
qubit coin space, He, spanned by two orthonormal vectors denoted {\R),\L)}. A generic 
state for the walker is 

oo 

l*>= Yl \x)®[a*\R)+b x \L)] (1) 

x=— oo 

in terms of complex coefficients satisfying the normalization condition Yl x \ a x\ + 1^1 2 = 1 
(in what follows, the summation limits are left implicit). 
A step of the walk is described by the unitary operator 

U = S ■ (Ip ® C) (2) 

where C is a suitable unitary operation in He and Ip is the identity in Tip. A conve- 
nient choice is a Hadamard operation, H\R) = (\R) + \L)) /y/2 and H\L) = (\R) - \L)) /y/2. 
When C = if, as in this work, one refers to the process as a Hadamard walk. The shift 
operator 

S = S R ® \R)(R\ + S L ®\L)(L\ (3) 

with Sr = J2 X \ x + an d $l = S R = J2 X \ x ~ l)( x \-> conditionally shifts the position one 
step to the right (left) for coin state R (L). Furthermore, it generates entanglement between 
the coin and position degrees of freedom. 

The evolution of an initial state |^(0)) is given by 

|tf(t)) = I7*|#(0)) (4) 

where the non-negative integer t counts the discrete time steps that have been 
taken. The probability distribution for finding the walker at site x at time t is 
P(x;t) = | (x\ty(t)) \ 2 = |<Xr| 2 + \b x \ 2 . The variance of this distribution increases quadrati- 
cally with time 9j as opposed to the classical random walk, in which the increase is only 



linear. This advantage in the spreading speed in the quantum case is directly related to 
quantum interference effects and is eventually lost in the presence of decoherence |8|. 

One of the early papers on quantum walks, due to Nayak and Vishwanath [10], has 
shown that Fourier analysis can be successfully used to obtain integral expressions for the 
amplitudes a x (t) and b x (t) for given initial conditions. The resulting quadratures can be 
evaluated in the long-time limit. The usefulness of this approach has been limited because 
the detailed calculations are lengthy and must be individually worked out for each initial 
condition. In this paper, we shall use the dual Fourier space to obtain information about 
the asymptotic entanglement, and so we review those results. 



Fourier transform 



The dual space Hk is spanned by the Fourier transformed kets \k) = XL e \ x ), where 
the wavenumber k is real and restricted to [— n, ir]. The state vector (CQ) can then be written 



dk 

— \k)® ~a k \R) + b k \L) 
An 



(5) 



where the k-amplitudes a k = (k,R\^f) and b k = (k,L\^f) are related to the position ampli- 
tudes by 

d k = ^e- ikx a x and b k = e~ ikx b x . (6) 

X X 

The shift operator, defined in (jHJ), is diagonal in 7i k space: S\k, R) = e~ lk \k, R) and 
S\k,L) = e lk \k,L). A step in the evolution may be expressed through the evolution op- 
erator U k in fc-space as 



\Q k (t + l)) = U k \Q k (t)) = -= 



I I e lk e %k 



!**(*)> 



(7) 



V% \ e ik —e i 

where \$ k ) = (k\^f) is the spinor (a kl b k ) T . This operator has eigenvectors {(p^ 1 ' 2 ^} given by 



Uk 



Vk 



Uk 



(8) 



where ak and /3 k are the real, positive functions 

1 

7! 



1 + cos 2 k — cos kVl + cos 2 k 



Pk = 



V2 



1 + cos k + cos kVl + cos 2 k 



-1/2 



-1/2 



(9) 



and 



u k = e ik 



Vk 



V2e~ iuJk - e~ ik (10) 



Wk = -^/2 e ^k_ e -ik 

The frequency Uk, defined by 

sin Jv 

smu k = ^=r- 0J k e [-7r/2, 7r/2], (11) 
V2 

determines the eigenvalues, ±e TtuJk , of Using the spectral decomposition for Uk, the time 
evolution of an initial spinor can be expressed as 

|$ fc (f)> = E/J|S*(0)> = (12) 

+(-i)V-^ 2) |$ fc (o)> bl 2) ). 

In principle, this expression can be transformed back to position space and the probabil- 
ity distribution P(x,t) can be obtained. This approach requires evaluation of complicated 
integrals which, for arbitrary times, can only be done numerically. However, in the long 
time limit, stationary phase methods can be used to approximate the resulting integrals for 
given initial conditions, an approach illustrated in [loj]. In this work, we bypass these tech- 
nical difficulties, because the asymptotic entanglement introduced by Uj., may be quantified 
directly from eq. (IT2"|) . without transforming back to position space. 



Entropy of entanglement 

Consider the density operator p = \^f)(^f\. Entanglement for pure states can be quan- 
tified by the von Neumann entropy of the reduced density operator p c = tr(p), where the 
partial trace is taken over position (or alternatively, wavenumber k). Note that, in general 
tr(p 2 c ) < 1, i.e. the reduced operator p c corresponds to a statistical mixture. The associated 
von Neumann entropy 

S E = -tr(p c log 2 p c ), (13) 

also known as entropy of entanglement, quantifies the quantum correlations present in the 
pure state p [ll| . It is zero for a product state and unity for a maximally entangled coin state. 
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It is also invariant under local unitary transformations, a usual requirement for entanglement 

measures HQ. 

The entropy of entanglement can be obtained after diagonalisation of p c . This operator, 
which acts in Tic, is represented by the Hermitean matrix 

Pc = , (14) 



where 

5 



5* C 



elk . . 2 

I^T = / o I 

-7T 271_ 



^ a x .&*. = y ^afc&fc (15) 



'ft I 



Normalization requires that tr(p) = 1 and A + C = 1. In terms of 

A = AC-\B\ 2 (16) 
the real, positive eigenvalues 7*1,2 of this operator are given by 



M 2 



1 ± Vl -4A 



(17) 



and the reduced entropy can be calculated from eq. (TT3l) as = — (ri log 2 7*1 + r 2 log 2 r 2 ). 
In the rest of this work, we shall be concerned with clarifying the asymptotic (long time) 
value of Se for both local and non-local initial conditions. 



Asymptotic entanglement from local initial conditions: a simple example 



As a simple application, consider the particular localized initial condition 

|*(0)> = |0) ® \L), (18) 

which implies eifc(O) = and 6fc(0) = 1. For this simple case, the spinor components at time 
t from eq. (fl2j) . have been explicitly calculated in Ref. 10] as 



a k (t) 
kit) 



/e 



ik 



2 Vl + cos 2 k 

1 / cos/c 

2 V + Vl + cos 2 k 



e 

cos k 



(19) 



Vl + cos 2 k 



Juj k t 



The relevant quantities for the entropy entanglement are A and B, defined in eqs. (ll5l) . 
After some manipulation, from ffl9l) we obtain 

™ dkl- (-l)*cos(2wjfct) 



A(t) 



2vr 2(1 + cos 2 k) 



r dk ie ik . , 

J-tt 27T 2V1 + cos 2 k 

, t f cos(2a>£;t) cos k . . . 
V v 1 + cos^ K 

The time dependence of these expressions vanishes in the long time limit and we obtain the 
exact asymptotic values, 

j , 1 f* dk 1 Vz 

A = lim A(t) = — = — (21) 

t-oo w 2 / , 2tt 1 + cos 2 k 4 v 7 



B = lim B(t) = - 



ff dk cos 2 A; .2 - v 7 ^ 



t-oo 2 J_ v 2tt 1 + cos 2 4 

From eq. ffl6l) we obtain 

A = A = ^-i (22) 

and the exact eigenvalues of the reduced density operator are r\ = l/V^andl-l/v 7 ^. These 
eigenvalues yield the asymptotic value for the entropy of entanglement So ~ 0.87243 . . ., 
in agreement with the numerical observations reported in Ref. [6|. When the procedure 
outlined above is repeated for arbitrary initial coin states, the asymptotic entanglement 
level is dependent on the initial coin state. Let us parametrize the initial coin in terms of 
two real angles a G [— n, tt] and (3 G [— n, tt], 

\x) = -4= {cosa\R) + e ip sina\L)) . (23) 
v2 

Then, the dependence of A, defined in eq (1161) . on the initial coin is 



A = A - 2b\ cos p sin(4a) , (24) 

with Ao given in eq. (122j) and b\ — (2 — v / 2)/4. We include the details of the calculation 
leading to this expression in Appendix A. Thus, A (a, /3) is a smooth, odd, periodic function 
of a and only in the special cases 
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FIG. 1: Panel (a): asymptotic entropy of entanglement, Se, for localized initial states as a function 
of a and /3, from eq. (|23p (color online). For /? = 0, (thick line), /3 = 7r/4, (thin line) and (3 = it/2 
(dashed line). Panel (b): time evolution of the entropy of entanglement, Se, for a = —ir/8 and 
(3 = 0, (full asymptotic entanglement), (3 = ir/2 (intermediate asymptotic entanglement) and 
(3 = 7r/4 (minimimum asymptotic entanglement). 

i) f3 — ±7r/2; for arbitrary a 

ii) a = 0, ±7r/4, ±7r/2, ±37r/4,±7r; for arbitrary /3 

the asymptotic entanglement is Se ~ 0.872. . .. In panel (a) of Fig. [T]the variation of the 
asymptotic entanglement associated to localized initial coins is shown. Full entanglement 
can be obtained for special initial coins. In panel (b) the time evolution of the entropy of 
entaglement is shown in three special cases with a = —3ir/8: for (3 = 0, full entanglement 
is achieved rather fast and the oscillations are supressed. For (3 = it/2 the intermediate 
asymptotic entanglement level 0.872 . . . results. For /3 = tt, the minimum asymptotic entan- 
glement is obtained, but the oscillations are larger and the convergence slow. If non-local 
initial conditions are considered, lower values for asymptotic entanglement may be obtained. 
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FIG. 2: Evolution of the entropy of entanglement for the initial delocalized sates from eq. (|25|) (up 
and down black triangles respectively) and for a localized state |0) ® \x) (gray circles). The time 
evolution was calculated from eqs. (|30p . The horizontal lines represent the asymptotic entanglement 
levels, 0.661, 0.872 and 0.979, obtained from eqs. (|3ip as explained in the text. The inset shows the 
first steps in detail. 

III. NON-LOCAL INITIAL CONDITIONS 

Most previous work on quantum walks has dealt with initial wavevectors localized in 
a position eigenstate |0). When non-local initial conditions are considered, new features 
emerge. Let us consider a quantum walk initialized in a simple uniform superposition of two 
position eigenstates such as, 

l*±) = ' " + ^ ® IX) (25) 

with the initial coin fixed at \x) = (\R) + i\L))/y/2. For this initial coin, any localized initial 
condition yields asymptotic entanglement Se ~ 0.872 . . .. The entanglement induced by the 
evolution operator when starting from these initial states is shown in Fig. [2j The asymp- 
totic values are S + ~ 0.979 and £_ ~ 0.661, respectively. Below, we provide an analytical 
explanation for these observed values. 

As can be seen in Fig. [21 the rate at which the asymptotic value for entanglement is 
approached is faster for higher asymptotic entanglement levels. The inset in this figure shows 
the first five steps in detail. The local initial condition j^o) is fully entangled (Se = 1) after 
the first time step and reaches its asymptotic level after ~ 10 steps. The non local conditions 
|\&±) have the same evolution for Se in the first two steps, but the phase difference causes 
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FIG. 3: Asymptotic entropy of entanglement, Se, as a function of A and \B\, denned in eqs. (115|) 

very different entanglement levels after the third time step: |^+) reaches its asymptotic level 
after three steps, while |^_) takes about 30 time steps to stabilize. 

Asymptotic entanglement 

Consider the problem of determining the asymptotic entanglement for non-local initial 
conditions of the form 

\*(0)) = '52c x (Q)\x)®\x) (26) 

X 

with real c x (0). The coin state \x) in eq. fl26l) is such that 

b x (Q) = ^(O) = X r or, equivalently, 6^(0) = ia k (0). (27) 
v2 

This restriction considerably simplifies the algebra, but it is not essential and the method 
applies equally well to arbitrary initial coin states. 

The eigenvalues of the reduced density operator depend on the real coefficients A, C and 
on the complex one B, defined in eqs. fll5l) . Fig. [3] shows the entropy of entanglement as 
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a function of these coefficients. Maximum entanglement would be obtained for A = 1/2 
and i? = 0, when the reduced operator corresponds to the minimum information mixture 
Pc = 1/2. 

In order to find expressions for A and B, we start by rewriting eq. f[T2l in the more 
explicit form, 

a k (t) = a\F k u k e-^ + (-1)'/3 2 GW^ 

~b k (t) = alF k v k e-^ + (-lYPlGkWke^. (28) 

Here, a k and fl k are the real, positive functions defined in eq. and u k , v k and w k , the other 
part of the eigenvectors of U k , are defined in ([TO]) . The dependence on the initial conditions 
is contained in the complex factors 

F k = u* k a k (0)+v* k b k (0) 

G k = u* k a k (0) + w* k b k (0)- (29) 

The required expressions for \a k \ 2 , \b k \ 2 and a k b k are 

\~a k (t)\ 2 = ai\F k \ 2 + Pt\G k \ 2 + (-lY2alPlRe[F k Gle- 2 ^ 

\b k (t)\ 2 = at\F k \ 2 \v k \ 2 + (3t\G k \ 2 \w k \ 2 + {-^alPlRelFkGlvkwle- 2 ^} (30) 
~a k (t)bl(t) = ai\F k \ 2 u k v* k + p A k \G k \ 2 u k wl + (-l) 4 ^ 2 [F fc G> fc ^ e - 2i ^ + F* k G k u k vle 2wkt ] . 

In the long-time limit, the contribution of the time-dependent terms in the k— integrals 
of eqs. (13"UI) vanishes as t _1//2 , as shown in detail in [h]]. The asymptotic values A, B and C 
can be obtained from the time-independent expressions, 




which hold for arbitrary initial conditions. 

Now we particularize these expressions for the initial states defined in eq. ( 1251) . Using the 
symmetry condition, eq. ( 1271) . the required squared moduli \F k \ 2 and \G k \ 2 can be expressed 
as 

\F k \ 2 = 4|fi fc (0)| 2 [l-cos(A;-cu fc + 7r/4)] 

\G k \ 2 = 4|5 fc (0)| 2 [1 + cos(A; + cu k + tt/4)] . (32) 
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Since A + C = 1, we need only perform the integration for A, 

A = fJ^ W = £§W)M°)| 2 (33) 
where Q(k) is a real function, 

Q(fc) = 4 [o£ (1 - cos(A; - u k + tt/4)) + $ (1 + cos(fc + u k + tt/4))] (34) 

which, after some manipulation, can be simply expressed as 

sin k cos 

QW = 1 + TwP (35) 

Note that this weight function satisfies Q(k) > and ^Q(k) = 1. 
The asymptotic form for I? is obtained from 

5 = fJ^*R= J[t m M0)|2 (36) 

whith given by, 

e~ lfc sin£; . . 

fl W = TwI' (37) 

The required quantity, Se, can now be evaluated from eqs. (l33]) and ( 1361) for different 
initial conditions satisfying eq. (12"7|) . 

For the local initial condition |0) ® we have 1 5^(0) | 2 = 1/2 and, from eq. (1531 . 

A = (7 = 1/2 results immediately. In this case, eq. (TPTj) reduces to fi j2 = | ± and 

the asymptotic eigenvalues are determined by \B\ alone. This quantity is obtained from 
eq. (l36l) . after noticing that if |afc(0)| 2 is an even function of k, only the imaginary part of 
R(k) contributes. We obtain, 

B =\B^ )\ = y^ (38) 



and the asymptotic eigenvalues, r\ = l/v2 and r 2 = 1 — l/y/2 give the asymptotic entan- 
glement level of So ~ 0.872. 

We now return to the case of non-local initial conditions |\P±) defined in eq. (I25p . for 
which 

cos 2 k for |^ + ) 

l a *(°) I = 2I £ (39) 
I sin k for 

Since J^. 7^Q{k) cos 2 A; = § Q(k) sin 2 jfe = §, we obtain from eq. ([33]), A = C = 1/2 

and the eigenvalues are also determined by \B\ alone. Since \dk(0) | 2 is still even, only 
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the imaginary part of R(k) makes a contribution. When inserted in eq. (1361) . these initial 
conditions result in 

B + = \B(* + )\ = \j: n iR(k) cos 2 £;| = (v ^~ 1)2 
R_ = |B(*_)| = | f£-R(&) sin 2 A; | ^-(v^-l) 2 . 

(40) 

Note that these values are related by B = | (£>_ + B + ). The exact eigenvalues are, 

r l = 2 - V2, r 2 = V2-l for |^+) (41) 

and 

ri = 2(^2-1), r 2 = 3-2^2 for |*_), (42) 
and the corresponding asymptotic entanglements are, 

S + = — - 1 - log^v 7 ^ - 1) = 0.97866 . . . (43) 
v2 

= -2(V2 - 1) (l + v 7 ^ \og 2 (V2 - 1)) = 0.66129 

respectively. These exact values are coincident with those obtained numerically, see Fig. [2J 
These are the maximum and minimum possible entanglement levels which can be obtained 
when starting in the position subspace TCx spanned by the kets | —1) and |1) with fixed initial 
coin. We show below that starting from a generic state in this subspace, all intermediate 
values of asymptotic entanglement are possible. 



Generic non-local initial state in TCi 

Consider a generic ket in Tii as the initial condition for position. We keep the same initial 
coin \x) which leads to a symmetric evolution in the local case. Thus, we consider initial 
states of the form 

\V(0,<p)) = (cosfl| - 1) +e-^sin0|l)) <8> |x>. (44) 

The parameters 9 G [— 7r/2,7r/2] and ip G [— 7r, 7r] are real angles. The initial amplitudes are 
5 fc (0) = (e ik cos6 + e ~ i{ - k+ ^ sinfl) /V2, b k (0) = ia k (0) so, 

\a k (0)\ 2 = M0)\ 2 = \ [1 + sin(20) cos(2A; + <p)\ . (45) 
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FIG. 4: Asymptotic entropy of entanglement Se(0, <p) as a function of the initial state from eq. (|44p . 
The entropy is calculated from the expressions for C and \B\ given in eqs. ()46p and (|47p . 



We use eq. (13311 to obtain, 

1 f n dk 

A = - + sin(2#) sh%) / — Q(jfe) sin(2A;) 

= --S + siii(20)siii(^). (46) 

For arbitrary y> and 8 = 0, ±7r/2 (indicating localized initial positions) or for arbitrary 9 
and ip = 0, ±7r (indicating relative phases zero or 7r between initial position eigenstates) we 
have A = C = 1/2 and the asymptotic entropy is determined by the value of \B\ alone. 

The expression for B obtained from eq. ( 1361) is, 

B = -B+ sin(2#) sin(^) - i{B - B' sin(2#) sin(^)) (47) 
with B defined in eq. (1381) and 



14 



<p 




FIG. 5: Contour plot the surface shown in Fig. [3] using the eigenvalues from eq. (|49p . Red areas 
indicate maxima and blue areas, the minimum values. 



The asymptotic eigenvalues r 12 are obtained from eq. (TT7j) as 

1 



r i,2(0,(p) 



± 



(B - B' sin(2#) cos(p)) + 2{B + sin(2fl) sin(p))' 



1/2 



(49) 



The resulting asymptotic entanglement Se{9, <p) is shown in Fig. HI A contour plot of this 
surface, Fig. [5J shows that there are only two maximum and two minimum points (for 
initial conditions |^+) and |^-) or (6, if) = (vr/4, 0) and (7r/4,7r)) for which the asymptotic 
entanglements are S + ~ 0.979 and S*_ ~ 0.661, respectively. 



More general non-local initial states 

In order to further illustrate the effects that the non-locality in the initial condition can 
have on the asymptotic entanglement level, we consider an initial Gaussian wave packet with 
a characteristic spread a 3> 1 in position space with the same coin state \x) as before. In this 
case, the Fourier transformed coefficients 5^(0) oc ^/cr e~ fc2 °" 2 / 2 correspond to a well-localized 
state in k space with lim^oo 1 5,^(0)1 2 = 2ix5(k), where 5(h) is Dirac's delta function. In this 
limit, the elements of the asymptotic reduced density matrix can be trivially evaluated from 
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FIG. 6: Asymptotic entropy of entanglement Se{0,^P = 0) for initial states in 7i\ (see eq. (|44p ). 
Notice that the 6 increases leftwise. This curve was generated using the eigenvalues given by 
eq. (|49p . with ip = 0. The maximum and minimum values (dotted lines) can be calculated with 
arbitrary precision from the exact eigenvalues given in eqs. (|41|) and eqs, (|42|) . The dashed line 
indicates the level of asymptotic entanglement associated to local initial conditions. 



the eigenvalues are r% — 1 and r 2 = and the corresponding asymptotic entropy of entan- 
glement vanishes. 

This simple example shows that for a particular uniform distribution in position space, a 
product state results in the long time limit. Of course, this is not true for arbitrary relative 
phases between the initial position eigenstates. However, it is clear that if more sites are 
initially occupied with appropriate relative phases, lower asymptotic entanglement levels 
may be obtained, until eventually a product state is reached. A more detailed analysis of 
this example shows that for cr 3> 1, the smaller eigenvalue approaches zero as (2cx)~ 4 and the 
asymptotic entropy of entanglement decays as Sg ~ log 2 <r/4cr 4 + 0(cr~ 8 ). Thus, for a > 10 
the asymptotic entanglement is already quite small, Se ~ 10~ 4 . Thus, if many sites are 
initially occupied, the entanglement level at long times becomes negligible. 



eqs. d33J and (JH]) 




16 



IV. CONCLUSIONS 



The long-time (asymptotic) entanglement properties of the Hadamard walk on the line 
are analytically investigated using the Fourier representation. The von Neumann entropy 
of the reduced density operator is used to quantify entanglement between the coin and 
position degrees of freedom. The fact that the evolution operator of a quantum walk is 
diagonal in fc-space allows us to obtain clean, exact expressions for the asymptotic entropy 
of entanglement, Se, for different classes of initial conditions. 

An expression for the exact value of asymptotic entanglement for localized initial states 
|0) <8> |x), with arbitrary coin \x) = (cos a\R) + e l/3 sin /a/2, has been obtained an- 
alytically. This entanglement changes between full entanglement and a minimum of 
Se = 0.736. . .. For full entanglement, the convergence to the asymptotic value is signifi- 
cantly faster than in the case of minimum entanglement. Lower asymptotic entanglement 
levels may be obtained if non-local initial conditions are considered. 

We considered in detail the case of initial conditions in the position subspace spanned by 
| ± 1), with fixed coin (\R) +i\L)/y/2, and obtained an exact expression for the asymptotic 
entropy of entanglement. This expression shows that it varies smoothly between the extreme 
values 5_ 0.661 and S + 0.979. As expected for this coin, the localized initial positions 
| ± 1) have an intermediate entanglement of So = 0.872 .... In order to explore the effect of 
increasing non-locality, the asymptotic entanglement the case of an initial gaussian profile 
in position space, with characteristic spread a ^> 1, is considered. For the particular phase 
relation considered, the resulting asymptotic entanglement decays fast with increasing initial 
non-locality. Thus, if many sites are initially occupied, a negligible entanglement level may 
be obtained at long times. 

The results presented in this work need to be extended to less simple systems. Most 
likely, quantum walks with either more particles, more dimensions or both will be required 
to be useful for algorithmic applications. The problem of entanglement in such systems is 
more involved. For example, in a quantum walk with two non-interacting particles there are 
four degrees of freedom and several kinds of entanglement may coexist. Some initial work 
in this direction is presently under way. 
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Appendix A: 

Asymptotic entanglement level for localized initial conditions. 



In this Appendix we obtain the analytical expression for the dependence of the asymp- 
totic entanglement level on the initial coin state for the case of initially localized position 
eigenstates. Let us consider as an initial condition a position eigenstate, which we take as 
as x = without loss of generality. Let us also consider a generic coin 

|$(0)) = |0) <g> 4= {cosa\R) + e ip sma\L)) (50) 
v2 

where a, (3 are two real angles. The Fourier-transformed initial coefficients, eq. ([6]), are 

5o(0) = a (0) = cos a 

6o(0) = 6o(0) = e i(3 sma (51) 
From expressions ( 129|) . for this class of initial conditions we obtain 

\F k \ 2 = cos 2 a + |ffc| 2 sin 2 a + sin(2a) Re [«^fce _;i/3 ] 

= F e (k) + F (k) 
\Gk\ 2 = cos 2 a + \wk\ 2 sin 2 a + sin(2a) Re [w^u^e - ^] . 

= G e (k) + G (k) (52) 

where the even and odd parts are the real functions 

F e (k) = cos 2 a + \vk\ 2 sin 2 a 

— (l — v / 2cos(u;fc — sin(2a) cos/3 
F (k) = — v / 2sin(u; fe — k) sin(2o;) sin/3 

G e (k) = cos 2 a + | w k | 2 sin 2 a (53) 

— ^1 + V / 2cos(co , fc + k) \ sin(2a) cos/3 
G Q (k) = — v2 sin(u; fc + k) sin(2o!) sin/3. 

The eigenvalues (TTTjl of the reduced density operator, which determine the entanglement 
level, depend only on its determinant 

A = AC- \B\ 2 = C- (C 2 + \B\ 2 ). (54) 
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The asymptotic entanglement will be independent of the initial coin state if A does not 
depend on a or f3. From the general asymptotic expressions (13T1) which are valid in the 
long-time limit for arbitrary initial conditions, we obtain 

C = c\ + (02 — ci) sin 2 a + C3 sin (2a) cos (3 

B = 61 + (6 2 — b\) sin 2 a + sin(2a) (63 cos (3 + i& 4 sin (3) 



where the real coefficients Cj ; bj ctrc explicitly 



V2 
4 



c 2 = 
h = c 3 



dk , a. 1 o ^4 1 i9\ \/2 
— (a fc v fc + A,K| ) = — 
_ Zix 4 



1 y/2 

2 ~ ~T 



r dfc 




L27 









- \/2cos(a; fc - fc)) + 

(l + v / 2cos(cu fc + fc) 

^1 — v / 2cos(cjfc — A;)j + 

/3fcK| 2 (l + >/2cos(wfc + k) 
4 (l-V2cos{uj k -k) 2 y + 

Pi (l + v / 2cos(cu fc + fc) 

2 / 2~7T H sill2 (^ - + 



1 

X ~~ 2 



1 v 7 ^ 

2 ~ ~T 



/^sin 2 (u; fc + A; 



V2-1 



(55) 



(56) 



These coefficients can be conveniently expressed in terms of one of them, i.e 61, since they 
satisfy the relations 

Ci = &i + §, c 2 -c 1 = -2fei, c 3 = 61, 
62 = -h, 63 = 61, 64 = \/2 61. 

Thus, eqs. ( )55l) can be expressed as 

C = bi 3 + v 7 ^- 2 sin 2 a + sin(2a) cos/5 

B = bi 1 - 2 sin 2 a + sin(2a) fcos/3 -H\/2sin/?J 

(57) 
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and the exact determinant is,, 



A = A - 26? cos sin(4a) (58) 

where 

A = c 1 (l-c 1 )-bj = y^ T ±. (59) 

These expressions can be used to calculate the exact asymptotic value for the entropy of 
entanglement for local initial states with arbitrary coins. 
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